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Abstract 

I analyze the lattice simulation data of the Domain Wall Fermion in quaternion basis. As 
pointed out by Atiyah and Ward, the minimum action solution for SU(2) Yang-Mills fields 
in Euclidean 4-space correspond, via Penrose twistor transform, to algebraic bundles on the 
complex projective 3-space. Assuming dominance of correlation between the fermions on the 
domain walls via exchange of instantons, I extract parameters necessary for defining gauge 
fields of Atiyah- Ward ansatz. The QCD effective coupling in the infrared and the relation 
between the number of flavors and the infrared fixed point is investigated. 

Consequence of this lepton flavor assignment to phenomenology of baryons is also discussed. 
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^The properties of infrared QCD depend on the flavor number of fermions. In 1982, Banks and 
^J2aks[T] showed that the 3-loop (3 function of an SU(3) gauge theory with Np Dirac fermions in the 
j^fepresentation R is 

/5o = II -^T{R)Nf 

(3i = 102- {20 + 4C2{R))T{R)Nf 

,285 5033,, 325 ,,2N 
132 = Np + Nl) (1) 

From the ratio of /3o and (3i which are invariant under the renormalization scheme, they showed 
that there appears an infrared fixed point when Np = = 16.5 

Using the staggered fermion and the Schrodinger functional scheme, Appelquist et al[21 E] 
claimed that if Nf is in the range of 12 < Nf < 16, the infrared behaviors are governed by the IR 
fixed point. In the Nf = 4 system, the step scalng function shows size dependence whe the coupling 
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constant u is large[3lll]. Since staggered fermion has the 4 times taste degeneracy, the interpretaion 
of the effective Nj is not simple. Recently, Fodor et al. [5] showed that the dependence of chiral 
condensates and the pion decay amplitude do not show the expected dependence on the QCD 
anomalous dimension 7 and cast doubt on the conformal window hypothesis. Staggered fermion do 
not have Dirac 7 matrices and the algebraic structure could be diffrent from that of Wilson fermion. 

Dietrich and Sannino[6] studied critical flavor number for a presence of Banks-Zaks type fixed 
point in two-loop /3 function for a generic non-Abelian gauge theory with fermionic matter in a given 
representation of SU(N). Non-perturbative corrections to these perturbative theory is not clear. 

Grunberg[71 El |9], assumed that there is two-phase structure in QCD as the number of flavors 
is changed. For < Nf < = 16.5, the conformal field amplitude D{Q'^) is calcuated from 
perturbation theory Dpt{Q^) and in the region Nf < Nj- the non-perturbative correction Dnp{Q'^) 
is to be added. It is also assumed that there is non-trivial negative UV fixed point auv The (3 
function which allows non-trivial auv and ajR is expressed as 

/3(a) = -/3oa' + M'----, (2) 

where $1 = pPq. The parameter p is adjusted such that (3{a) = for /3o = i.e. Nj = 16.5 and for 
each Nf 

/3(a) = -/3oa^ + (3) 
coincides with the true /3{a). He obtained am^uv as 

am,uv = ±~e{l + l^~e + ---) (4) 
^ P20 

where e = {~'^y^'^ and /3jo are the /3j at Nf = 16.5. In this modified Banks-Zaks fixed point 
theory, the conformal window is predicted in 4 < Nf < 6 
The resummed running coupling Ro{Q'^) is defined as 

Ro{Q') = l^a{e)<p{e/Q')^ (5) 

where a{k'^) = ""^^^ ^ and is observable dependent Feynman integrand. 

Singulaity of the Borel transform of Rq defined as 

POO 

R^{Q^) = / B{z)e-'/^^'^'^dz (6) 



~'0 

on the real axis in the z complex plane are called renormalons. One assumes 

where 5 = ^ and z„ = pi (5q is the renormalon location. Banks-Zaks expansion with infrared fixed 

Po 

points has factorial divergence. 

Although it is not a direct evidence, infrared freezing of the QCD effective coupling is an in- 
dication of proximity of the system to the conformal window. In MOM scheme, the quark-gluon 
coupling as{q) in the Coulomb gauge of the 2-1-1 flavor domain wall fermion (DWF)[T2], and the 
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Figure 1: The effective coupling of QCD. Jlab 
experiment and lattice simulation. Small disks 
from logioQ from 0.5 to 0.3 and the solid 
curve are the <ygi{q) of JLab. Diamonds and 
squares are as{q) of ghost-gluon coupling of 
quark mass = 0.02/a and m„ = 0.01/a, 
respectively [TT]. Three large blue disks are 
as{q) of the quark-gluon coupling. 
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Figure 2: The beta function obtained from the 
fitted function in the Figure 1. 



ghost-gluon coupling as(q')[TT] are calculated. In Fig.l, the results are compared with the exper- 
imental extraction of C(gi{q) derived from the polarized electron proton scattering by the JLab 
group [13] • 

The fitted line of the (ygi{q) is derived from [13] and the corresponding beta function is calculated 
numerically using Mathematica[35j . which is shown in Figj2j In [18j, af'^^{q) of AdS^xS^ [19j theory 
at high energy and the phenomenological (ygi{q) at low energy are combined. In the infrared, the 
correction from af^"^^ is small. The coupling constant in different scheme for calculating a physical 
quantities are related by commensurate scale relation[15]. It is remarkable that oi^'^^^q) is close to 
agi{q). 

The beta function does not have zero, probably due to the isometry of AdS space mapping, 
— )• and the holographic coordinate space mapping z ^ Xz. 

The freezing of the running coupling is observed also in ar{s) in the infrared[T7]. In this work, 
essentially only the experimental data of r decay was used. Beneke and Braun[TB] pointed out 
crucial roles played by renormalons in the Brodsky-Repage-Mackenzie(BLM) prescription of scale 
fixing. In BLM the scale changing to Q* by incorporation of fermion loops has the same effect 
as resummation of diverging radiative corrections. The infrared freezing of the running coupling 
observed by Brodsky et al could be understood in the same framework as Grunberg's effective 
charge theory [TO] . 

Brodsky and Shrock[Tl| discuss low energy QCD phenomenology and point out that instanton 
effect may be suppressed by a presence of confinement induced gluon momentum kmin 

Lattice simulations of the effective QCD coupling in the MOM scheme, in contrast to that 
in SF scheme, suggest that the three flavor system is close to the conformal window. In the SF 
scheme [22]. the winding number of the system is incorporated through the Chern-Simons actions 
on the boundary field . It is, however, not clear how the topological winding number of different 
scale steps are matched . 

In the DWF, I define the phase of the quark wave function in the 5th dimension, and adjust 
the phase so that the wave function on the left wall and that on the right wall are correlated by 
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self-dual instantons. Here, the eigenstate of the triality is considered, but not the winding number 
of the Chern-Simons action. I guess three times more flavor configurations as compared to that of 
MOM scheme were required in the SF scheme. 

In order to justify this conjecture, I first discuss the instanton or self-dual gauge field and their 
treatment in quaternion basis. Importance of an instanton in the low energy QCD was pointed 
out by t'Hooft in 1976 123], and Atiyah and Ward[25j showed that it emerges from the minimum 
action solution for SU(2) Yang-Mills theory. They applied the Penrose twistor approach |26j to the 
instanton problem and showed that there is one to one correspondence between a) self-dual solutions 
of the SU(2) Yang-Mills equation on 4-sphere up to gauge equivalence and b) isomorphism classes 
of 2-dimensional algebraic vector bundles E over projective 3-space P3 satisfying two conditions: 

1) E has a symplectic structure and 

2) The restriction to E of every real line of P3 is algebraically trivial. 

The DWF contains the 5th dimension which could play the same role as z in AdS/QCD. It 
preserves the chiral symmetry in the mass- less limit. We performed a simulation using the DWF 
gauge configuration of the RBC/UKQCD collaboration by fixing to the Landau gauge and to 
the Coulomb gauge. 

In [21], I analyzed the DWF simulation data using quaternions. The SU(2) spin algebra is 
expressible by quaternions H invented by Hamilton, which is a generalization of the complex number 
q = w + ix + + kz, endowed with the multiplication rule 

i2=f = k2 = -l, ij = k=-ji, jk = i = -kj, ki=j = -ik (8) 

The quaternion q = w + ix + ]y + kz can be represented by complex 2x2 matrix 

, f w + iz ix + y \ 
w + ix + w + kz = { . (9) 

y tx — y w — iz J 

In three dimensional space and in seven dimensional space R'', the cross product of two basis 
vectors is expressed by another basis vector. In the case of R'', the rule is given by 

ei X e^+i = ei+3 (10) 

where indices are permuted cyclically and translated modulo 7. The cross product in R^ is invariant 
under rotation of SO (3), but the cross product in R^ is not under rotation of 5*0(7), but under G2 
which is a sub-group of S0{7). Just as a, b, c in R'^ defined the quaternion algebra H, a, b in R^ 
and the cross product defines an octonion algebra O. 

The content of this paper is as follows. In sect. 2, I review the DWF propagator. I explain, in 
sect. 3 the Cartan's theory of spinor. The theory of Atiyah- Ward ansatz for the quark propagator is 
shown in sect. 4 and its application in the DWF lattice simulation in MOM scheme is summarized 
in sect. 5. The implication of triality in lepton flavor phenomenology is given in sect. 6 and on the 
B-meson leptonic decay etc. are given in sect7. 

2 DWF lattice simulation 

In the DWF method[l5], the mass of the fermion is 

m{s) = M sign{s) (11) 
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The continuum fermion propagator is constructed from 

Dp = iltid^i + i'jbds + im{s) (12) 

When a plane wave is assumed in the first four cordinates, the effective one dimensional problem 
becomes 

Dp = h^d^ + im{s) + 7 • p (13) 

It is a linear ordinary differential equation of first order and the solution ip = U{s, X,p)u{X,p,a) 
where U{s, X,p) is a 4x4 matrix and u{\,p, a) is an s independent spinor, becomes 

^ / cosh{sK+) + k;S5[^(7 ■p-X)-m] sinh(sK+) s > 
^ ' \ cosh(s«;_) + «;I^75[i(7 ■ p — X) + m] sinh(s/«_) s < ^ ' 

where k± — {iX ± M)^ 

The eigenvalue equation on the lattice in the s space becomes 

^^,[Uis + 1) - U{s - 1)] + ^[U{s + 1) + U{s - 1) - 2U{s)] 

+ (^Ylln sm{pi,) - F{p) + m{s) + iA j U{s) = Q (15) 

where F{p) = E^=i[l " cos(p^)]. 

The positive chirality component U^{s) and the negative chirality component U-{s) satisfy the 
relations 

4 

U+{s + 1) = [1 + F{p) -iX- m{s)]U+{s) - I J2 a^sin(p^)^7_(s) (16) 

4 

C/_(s - 1) = [1 + Fip) -iX- m{s)]U4s) - i J](7^sin(p^)C/_(s) (17) 

U{s + 1) = r+C/(s), s>0 (18) 
U(s-l) = T~U{s), s<0 (19) 

rj.+ _( 1 + ^(P) - - M -i a, sin(p^) \ 

~ a sinfo 1 ^-^^ ^ ^ 

where G(p) = E^=iSin^(p^) 

= 4^?^'^-»'(s, s';p) = (21) 

2 7 ■ p 

This expression shows infrared singularity. This singularity comes from the bound state pole in 
the lefthanded quark propagator. 

GAp) = ^ , Gnip) = ^ , (22) 
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Figure 3: Left handed chiral fermion amplitude 
A/, (p) (blue) and y4i(p)p^ (green). M = 1 
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Figure 5: Negative mass chiral fermion ampli- 
tude Aji{p) and Ar{p)p'^ {green). M = — 1.8 

When M > and s > > s' 
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Figure 4: Right handed chiral fermion ampli- 
tude Ar{p). M = 1 
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Figure 6: Negative mass chiral fermion ampli- 
tude v4L(p)(blue). M = -1.8 



[Gl{p)]ss' = Ai(p)e-+(P)^+--(P)«' (23) 

has the singularity, which is however cancelled by the gauge part [IHl HZ] ■ 
The kinetic energy extracted from the propagator at s = s' = becomes 

fi — 1 

where Al{p)p'^ can be interpreted as the normalization factor of the zero mode. It is finite at p = 
as shown in FigJHl Here I took M = 1. 

When the domain wall height M is taken to be negative M = —1.8 as assumed in[T2], the 
amplitude AR{p)hecomes singular in the infrared. 



3 Cartan's theory of spinors 

The quaternion space spanned by H spanned by k} can be extended by introducing a new 

imaginary unit / and make the octonion spsce O = H + H/ spanned by {1, k, I, il, ij, ik}. 

Cartan studied 2i> dimensional Euclidean space E21, and semi-spinor C,a with components of odd 
or even number of indices are zero. He defined as those of even number of indices are and ip as 
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those of odd number of indices are 0, and arbitrary //—vector X(^,y). Triality transformation rotates 
24 dimensional bases defined by Cartan[ 



{^0, ^1, ^2, ^3, ^4}, {^12, ^31, ^23, ^14, ^24, ^34}, {^123, ^124, ^314, ^234, ^1234} 

{/v»l /v»2 3 /v»4\ I t'^ t*^ -t*^ 1- (^^\ 

The trihnear form in these bases is 

F = (pCXilj = (^12^14 - 61624 - 64623 + 62346) 
(63624 — 62634 — 64623 + 62346) 
+3^^(61634 — 63614 — 64623 + 62346) 

"64634 — 64614 — 64624 + 62346) 
+^'^ (—6634 + 636 ~ 646 + 646) 
+a;^'(-66i4 + 616 - 646 + 646) 
+a;^'(-6624 + 626 - 646 + 646) 

+x'' (6623 - 636 - 616 - 626) (26) 

There are three semi-spinors which have a quadratic form which is invariant with respect to the 
group of rotation 

$ = *0C</) = 66234 — 6364 — 6164 — 6264 

^ = VCV- = -6634 - 6614 - 6624 + 6623 (27) 

and the vector 

F = x^x^' + x^x^' + x^x^' + x^x^' (28) 

The trihnear form J-" is invariant and the three quadratic forms interchange, under transforma- 
tions which are classified as G(23), ^(12) and G(i3), which are shown in the Appendix . The product 
*(G(i2))G(i3) generates G(i32) and *(G(i3))G(i2) generates G(i23)- 

The semi-spinor 6? 6i) 6im and 6234 are devided into two classes. The first type satisfies 

^ = 66234 — 6264 — 6364 — 6164 = (29) 

and the second type satisfies 

^ = -6634 - 6634 + 6624 + 6623 = (30) 
With a use of quaternion bases l,i,j,k, the spinors and Ccp = ip' are defined as 

A. c \ c ^ \ c A \ c u I 6 + ^64, ^64 — 64 A /on 
= 6 + 64* + 64J +64^ = _L^- ^- ,r (^31) 

= 6234 - 63* - 6iJ - 62fc = ( ^'l' - 5"' 7^''^S'' ) (32) 

\ —K23 — «1 U234 + ^'512 / 



^The definition in [53], G(i2)G(i3) = G(i32) and G(i3)G(i2) = G(i23) are misprints. 
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Similarly, and Cip = ip' are defined as 



When two spinors ip and C'?/' or and C0 are combined together, they are expressed by octonion, 
and the automorphism in the space of octonion is G2 group, which is exceptional Lie group with 
14 dimensional representation. In G2, there is a specific automorphism, which is called triality. In 
the following subsection, I study the structure of the G2 group. 

The trilinear form ^(pCXip which is invariant under triality transformation contains the term 
2^^64623, 3:^64^123 and x^^34^i23 which transforms by applying das to -634a;^'^o, -6142;^' ^0 and 
—^24^:^'^ Co, respectively. 

The vector and are Pliicker coordinates of the line which joins the two points and ip' of 
quark specified by and quark specified by ip'. 



4 The Atiyah-Ward ansatz in the quark propagator 



The null twister ^ defined by Penrose is a combination of two spinors A = (Aa) and fi = (yu") in the 
complex projective space CP'^ space. It satisfies 



/ 1 \ 

1 

10 

\ 1 / 

Aa/i" + fi^'Xa = 



A*, A*) 



V Ai y 



The components of the null twister satisfies 



(35) 



(36) 



An open subset of the twister space PT defined as PT' = {[^] G PT|^ = (/i. A), A = (Aq, A^) 7^ 
(0, 0)} allows a mapping on E 



y -z 

z y 



/i" -A 
A^ AO 



(37) 



The Atiyah-Ward ansatz for the quark propagator was adopted by Corrigan and Goddard[ 
in the transformation function g{uj,Ti) where u and vr are complex two-spinor which satisfy in the 
CP^ space, 

5f(Aw, Avr) = (7(cj, vr), detg = l. 
When X = — ix ■ (7 is a quaternion and u = xtt, g{u, it) can be expressed as 



g{xTT,TT) = h{x,C)Hx,C) ^ 
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where ( = — , h{x, () is regular in |C| > 1 — e and k{x, () is regular in |(^| < 1 + e. 
In [37], the Ansatz 

^° = ( ""o e° ) ( ) ( ) 

- ( "o^' i^^ ) <-) 

was proposed as the transformation matrix. Here j = fi — u. 
In [31], a solution of equation "^ooQ = ^o? where 

^0 = f M , ^oo = f ] , (39) 

V '^ai y V v^2i v^22 y ^ ^ 

under the boundary condition \l/(a;,oo) = 1 was considered. In our system, I identify \l/o and \l/oo 
as gauge configuration on the left domain wall and the right domain wall. 



5 The quark propagator in Coulomb gauge and the QCD 
effective coupHng 

I first perform a gauge fixing in Landau gauge using the over relaxation method and then fix to the 
Coulomb gauge without touching the A4 component [20]. I then measure the quark propagator and 
fix the gauge in the space of the 5th dimension between the domain walls. 

Since we adopt the representation of fermion spinors such that the 75 is diagonal, the fermion 
field ip can be gauge transformed in the fifth dimension as 

^ e'"^^'^, ^ ^e-^''^^ (40) 

such that on the left and on the right domain wall, the propagator becomes approximately real by 
choosing pTj 

be minimum. 

The equation has two sets of solutions, which can be obtained numerically. I assign 

^o,L/n =(1' and v,,^_,,/, = ( ) (42) 

V Co do y^/^ V CL.-i dL,-i 

from the color diagonal components of the lattice data. 

Using the lattice wave function ip{p, k) and the plane wave x(P; k), where k is the coordinate in 
the 5th dimension, I define 

Tr{xip^, 0)^L/RiPx, 0))o-i + Tr{xiPy, 0)^/'l/r(Pj/, 0))cr2 + 
Tr{xiPz, 0)^pL/RiPz, 0))a-3 + Tr{x{pu ^)'^l/r{.Pu 0))«/ 

oc a°(Ti + a°(T2 + a^cTa + aP^il (43) 



and 

Tr{xiPx, Ls - l)i>L/RiPx, Ls - l))ai + Tr{x{py, - VjipL/niPy, ^ - l))cr2 + 
Tr{xiPz, Ls - 1)iI)l/r{Pz, ^ - l))a-3 + Tr{xipt, Ls - 1)iPl/r{,Pu Ls - 



oc a. 



cTi + a^^^Va + a^^Vg + a^~^i/ (44) 

If the configulation on the left domain wall and the right domain wall are correlated by self dual 
gauge field, the transformation matrix g{'~^, C,) gives a relation, 

ais-i bi.-i \ f Ce'^ f \ _ f Ce-'^ f \ f ao bo 



CL^.i dL^.i I V C'e-' J V C^'e^ A Co 



(45) 



where/ =/(,.-!). 

In our 5-dimesional domain wall fermion case, the introduced phase in the 5th direction ir] is 
incorporated in n and —irj is incorporated in u. 
The function / and / taken in ^37j is 

dne'^ — ■:r-^ — e"^ — ar,e~^ 

/ = f-^^ ■ («) 

In general Cq and Ci^_i are polynomials of (. I define 

^ = c_iC~' + c,C' + 6 (47) 
where c_i = cl^-i and ci = cq and 5 is a constant. They are defined by minimizing the difference 

a-Ls-i bi.-i \ { C^e'^ f \ f C^e"'^ / \ f a-o bo 



^^/^nc; oc-c-^j -To rVj(c: <«> 

It is the condition for the left wall and the right wall are correlated by the self-dual gauge transfor- 
mation. 

Since C is a parameter in CP'^ space, I adopt an ansatz (^e'^ = \J ^ ~^d^^^^^ derived e*^ 
and the 6 by solving the simultaneous equation using Mathematica|35] 

aL^.J + 6z..-iC-'e-^ = boC'e-'^ + dof 
CL,-if + dL^.iC^e-^^ = doC^e'^ 

From the numerical practice, I observe that the deviation /\L/R becomes small when \5\ is large, 
i.e. when ~ 5 ~ y is large. When 5 is large 

^ = c_.c-' + .^ + c.c = \ - <^-'^" ^ ^'^*J ^ i'^-'^" ^ ^ ■ ■ ■ 

and the Laurent series of / converges fast. 

In the form factor calculation I use the overlap of gauge transformed x and the plane wave at 
the middle of the domain walls. 

Corresponding to the sample-wise selection of larger absolute value of 5, i.e. a large zero-mode 
component, I assign a parameter ind = 1 for larger and 2 for smaller, and multiply the phase (— l)*"-*^ 
to the expectation value of the quark wave function that is used in the calculation of the effective 
mass of the quark. This phase factor can be absorbed in the phase e^^ introduced to approximately 
satisfy the U(l) real condition. 
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6 Lepton flavor phenomenology 



The mass of a fermion is produced by Yukawa term which couples the boundary layer s = iV of a 
quark of a chilarity and the boundary layer s = of a quark of an opposite chilarity. In the minimal 
supersymmetric standard model(MSSM), the left-handed quark is defined as 

- ( I:: ) (») 

which have SU{3)c x SU{2)l x U{1)y quantum numbers 3, 2, 1/3. The right-handed quarks ur and 
dpi are SU{2)i singlet and they have quantum numberd 3, 2,4/3 and 3, 2, —2/3, respectively. 
In the MSSM, the Higgs field is defined as 

H={%^ (51) 



and 



- ( 1?;'; ) - ( I. ) (52) 



The vacuum expectation value of H is 



{H^H) = (53) 
The down quark mass is produced by the Yukawa coupling 

y^^Q^E^i^ + h.c. (54) 

and the up quark mass is produced by 

yu'^Q^iT^n^^'^un + h.c. (55) 
The masses of m, c, t quarks m„ c,i and c?, s, h quarks f, are expressed as 

"^«,c,t = VuDu.c^u rnd,s,b = ^dyd,s,b (56) 

When the spinors and ip are assigned to a lepton, e, /i or r, the electro-magnetic interaction 
requires they belong to the same family G(23), G(i2) or G(i3). The quark interaction with neutrino 
is expected to be blind to the lepton flavor, and the neutrino mixing occurs. 

Leptons are characterized by the flavor triality. Non-Abelian discrete flavor symmetry is 
assumed to be broken to Z3 in the quarks and by assigning to = e^'^*/^ the Left-handed quarks and 
leptons are classified [301 EI] 

M, c?, e ~ 1 c, s, ~ t, 6, r 



c, s, ~ w t, 0, r ~ w 



Ue Ur Ug Ub \ f 1^/, Or Cg Cb \ f Vr tg U 

e~ dr dg db ) ^ \ fi~ Sr Sg Sb ) ^ \ hr bg bb 



and the Right-handed quarks and leptons 



u'^ , d'^ , e'^ 1 c"^ , s'^ , fi'^ ^ u , b'^ , t'^ oj"^ . 



e Ur Ug Ub \ ( Cr Cg Cb \ f T tr tg tb 

R \ Sr Sg Sb ) p, \ br bg bb y ^ 
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There are self-dual gauge vector fields in the color-isospin space and in color-spin space. The 
vector field has the SU{5) D SU{3)c x SU{2)iy x U{1) symmetry. In the standard model, W, Z 
and electromagnetic field A*^™ are unified in Weinberg-Salam model. The gluon fields are denoted 
Ai, A2, A3 and Aq. 

r w+ w~ z \ f w+ w- z \ f w+ z \ 

V Ai A2 A3 Aq \ A2 A3 Aq \ Ai A2 A3 Aq ) 

In the following figures, I assign the self dual gauge fields or A as 

(rf^ /-yi^ /y^ \ 

■^1 -^2 -^3 -^4 1 
Xi X2 X3 X4 J 

which has their specific triality. 

The quarks q = u,d, s, b, c, t has G2 3 SU{3)c x SU{2)spin symmetry, but triality blind. 

Since vector fields are defined from spinors as Pliicker coordinates, the gluons that interact with 
quarks but not with e, /i and r, the triality symmetry for gauge bosons A^"^, W and Z and gluons 
could be qualitatively different. In the case of a color-isospin triality, electron number, muon number 
and tau meson number are well defined, but u, d, s, c, b and t quarks mix with each other. I expect 
gluon-weak current coupling preserves flavor triality and color-spin triality, but the quark-gluon 
coupling is blind to flavor triality [30] and to color-spin triality. The absence of flavor changing 
neutral current could be due to the triality selection, although the GIM mechanism[33] could also 
play a role. 

The lattice simulation is the most efficient method to clarify the mechanism of confinement and 
chiral symmetry breaking, since they are essentially non-perturbative phenomena. 

In astrophysics, most astronomical data are taken from electro-magnetic detectors, and electro- 
magnetic wave from a triality family that is different from that of the electro-magnetic probe would 
evade the detection. 

7 Discussion 

I studied possible roles of triality in color-spin space and its consequence in the IR-QCD, and 
extension of flavor-preserving self-dual gauge field to flavor changing gauge boson. 

It is not evident that the bispinor system proposed by Cartan actually explains the lattice 
simulation results. However, the difference of a critical flavor number in SF scheme, which is three 
times larger than that in the MOM scheme could be understood, if the boundary condition imposed 
in SF scheme, in which at each recursive renormalization step, matching large L lattice data to small 
L lattice data is performed, requires more freedom in the functional space of fermions. 

From the symmetry of octonions, Cartan constructed the vector flelds xi,X2,X3,X4 and their 
duals x[, X2, x'^, x'^ from spinors. Since I consider a self-dual vector fleld that couple to spinors, I use 
X\^X2^X3 for the self-dual gluon flelds or flelds. From the invariant trilinear form, I consider 
the Coulomb potential and the transverse component. I omit diagrams in which X4 is exchanged, 
since in magneto static QCD|12], one can choose Aq = 0, and ignore x^ exchange in self-energy 
diagrams. 

The three loop self energy could play an important role in magnetic mass problem[l3] in flnite 
temperature QCD. 
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gram: g644a 



Figure 10: Coulomb potential self energy dia- 
gram: g644b 
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Figure 11: Coulomb potential self energy diagram: g644c 
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Figure 12: The decay of 5 to r u. Transverse Figure 13: Transverse polarization tree dia- 
polarization tree diagram la. gram 2a. 
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Figure 14: Transverse polarization tree dia- 
gram lb. 



Figure 15: Transverse polarization tree dia- 
gram 2b. 



The model of vertex correction via self-dual gauge fields can be applied to a modification of the 
SM in weak decay process B~ — )■ r~i/ |39], and other meson decays. Some clear deviation from the 
SM is reported in Penguin diagram dominated pro cesses f38]. In Figs JT2|T3|14|15[ I show the B~ 
decay into t~u in a direct process and in Figs JT6l[T7|18|19l in Penguin-like process. Penguin of SM 



respectively. The self-dual gluon exchange contains infrared divergence and non-perturbative sim- 
ulation is necessary. 

Georgi[27] assumed that a world with "unparticle" modifies the structure of the infrared fixed 
point. There are investigation of supersymmetry inspired QCD beta function with unparticle 
incorporated If the electron or muon in the detector selects one triality in their electro- weak 
interactions, most results of SM are not affected. The q — q pair in different triality sector could 
behave as the unparticle. 

In order to explain the mass of W and Z bosons and other fermions, dynamical electro-weak 
symmetry breaking was proposed by several authors [IHl E] with an introduction of technicolor. 
Although the calculation is not done, I guess the flavor changing self-dual gauge field exchange, i.e. 
xi,X2 or xs in Fig.3 and 4 are identified as self-dual gluon or W^, affects the mass of gauge bosons. 




and 





Figure 16: The decay of 5 to r u. The pen- 
guin diagram la. 



Figure 17: The penguin diagram 2a. 
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Figure 19: The penguin diagram 2b. 



It is necessary to extend the lattice simulation to larger lattice and study the continuum limit. 
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Appendix: Matrix elements of the triality transformation of the g2 algebra. 

In the table 1 indicates -1. In Gu, I replaced ^4 — )■ —^1234 to ^4 — )■ —^123 written in [23] (probably 
typos). 6*136*12 does not agree with G(i23) of [23]. We think ^4 — )■ — there should be replaced by 
^4 — )■ — and we define ^{G^^J^is = ^123 and *(G'^3)G'i2 = 6*132, which correspond to G(i23) and 
G(i32) respectively of [23] . 

^723 ~ ^132 = ~25gQj.^ and Gi^2 ~ ^123 = 25x^go' i-^- inversion is same as interchange except the 
sign of ^0 x'^ transformation. 

Under these triality automorphism, the set {F, $, \E'} is invariant. 
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